A single ion as a shot noise limited magnetic field gradient probe 
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It is expected that ion trap quantum computing can be made scalable through protocols that 
make use of transport of ion qubits between sub-regions within the ion trap. In this scenario, any 
magnetic field inhomogeneity the ion experiences during the transport, may lead to dephasing and 
loss of fidelity. Here we demonstrate how to measure, and compensate for, magnetic field gradients 
inside a segmented ion trap, by transporting a single ion over variable distances. We attain a 
relative magnetic field sensitivity of AB/Bo ~ 5 ■ 1CP 7 over a test distance of 140 u.m, which can be 
extended to the mm range, still with sub [im resolution. A fast experimental sequence is presented, 
facilitating its use as a magnetic field gradient calibration routine, and it is demonstrated that the 
main limitation is the quantum shot noise. 



PACS numbers: 42.50.Dv; 03.67.-a; 37.10.Ty 



I. INTRODUCTION 

Over the past decade, trapped ions have emerged as 
one of the most promising systems to realize quantum 
computation, with many key ingredients implemented ex- 
perimentally, such as single and multiple qubit gates [l[. 
It appears difficult however, to scale single zone traps 
significantly beyond the order of ten qubits Q. One 
way to store and manipulate such complex systems on 
a larger scale with a larger number of ions, is to divide 
the ion trap quantum computer into several regions - 
each holding only a few ions - combined with transport 
of the ions between them (3|. In order to realize this, 
it is paramount to make the transport between storage 
regions and processor work with minimum loss of coher- 
ence, which has proven to be a challenge so far. In this 
paper we investigate the transport of single ion qubits 
and address one of the leading dephasing mechanisms, 
arising from magnetic field gradients. Previous transport 
demonstrations [J] have avoided magnetic field induced 
dephasing by adding an extra set of qubit levels, which 
are less sensitive to field changes during transport. This, 
however, adds extra complexity to the scheme, and also 
may not work for all ion species. It has also been shown 
that adding extra identical transport sequences after a tt- 
inversion can cancel out unwanted phase shifts Q even 
when transporting across junctions, however, such extra 
phase-changing operations does not appear to be feasible 
when transporting qubits that are a part of larger en- 
tangled states. Another approach is to use two physical 
ions in an entangled state to form one logical qubit in a 
subspace which is free of decoherence from the magnetic 
field variations over the transported distance [f|, at the 
expense of needing twice as many ions. Even in this sce- 
nario however, the difference in magnetic field magnitude 
between the two paired ions needs to be sufficiently low. 

Here we present a technique for measuring, and com- 
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FIG. 1. Part a) shows a reduced level diagram of the 40 Ca+ 
ion with the levels that are relevant for the presented experi- 
ments, b) shows the spin echo pulse sequence used to measure 
the magnetic field gradient. The preparation and the readout 
each consists of several steps, as described in the text, and 
the dashed pulse is the electric potential change that causes 
the transport, while the solid boxes are laser pulses. 



pensating for, the magnetic field gradient inside the trap, 
using a single ion as a movable probe. We carefully 
tune the currents of external Helmholtz coils for very 
low B-ficld gradients. Stationary two ion crystals also 
allow for local measurements of gradients Q, although 
the transport method presented here adds further ad- 
vantages: The attainable accuracy of the measurement 
is higher, due to the fact that the ion can utilize a much 
larger distance over which changes to the field are picked 
up. The gradient is measurable over a much larger spa- 
tial extent with a sub-u.ni resolution, and extends to re- 
gions outside those accessible by the laser beams and the 
ion observation, for example into storage regions or even 
across corners and junctions. 

For general scientific interest we envisage applications 
of our method where the combination of high spatial res- 
olution and shot noise limited performance is used to 
map out the magnetic field even of single protons Q or 
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molecules @. We note that urn resolutions have also 
been reached using Bose-Einstein condensates [loj . 



II. EXPERIMENTAL SETUP AND METHOD 



In our experiment, we employ a segmented micro- 



structured Paul trap [111 ], where a combination of static 
and dynamical electric fields provide confinement for a 
single 40 Ca + ion, and where the segmented design of the 
trap allows for transport [l2[ • The strengths of the trap- 
ping potentials lead to harmonic oscillation frequencies 
of uj/2-k = {1.35,2.4,3} MHz in axial and both radial 
directions, respectively. The relevant internal states of 
the ion are shown in Fig. QJ,, where Dopplcr cooling is 
performed on the S1/2 — > P1/2 (cycling) transition near 
397 nm. The two ground state spin levels, \m,j = ±1/2), 
are used as qubit states, and operated on by employing 
stimulated Raman transitions using a separate 397 nm 
laser, which is detuned by A 2ir ■ 40 GHz from the cy- 
cling transition. In addition, a laser at 729 nm provides 
state selective shelving during readout, while lasers at 866 
and 854 nm prevent trapping in metastable states. Two 
coils outside of the vacuum chamber in a Hclmholtz-likc 
configuration are set up to give a B-field axis at a 45° an- 
gle with respect to the trap axis. During the experiment, 
a magnetic field offset is kept around ~ 0.7 mT giving 
a Zeeman splitting of the spin states of about 18 MHz, 
while gradients are tuned by the application of different 
currents to the two coils. 

For determination of magnetic gradients we employ the 
following sequence of operations, depicted in Fig. [T|d: i) 
Preparation pulses consisting of Dopplcr cooling pulses 
as well as optical pumping pulses initialize the ion to the 
rxij = +1/2 state, ii) spin echo sequence, with the trans- 
port pulse inserted between the first 7r/2-pulse and the 
7r-pulse. iii) readout pulses, which are carried out by an 
efficient state selective transfer on the 729 nm transition 
followed by fluorescence from the Doppler laser [l3[ 

Precise knowledge of the transported distance as a 
function of time comes from a calculation of the potential 
minimum of the trap for each value of an applied volt- 
age by means of a numerical field solver [l4j]- In ref. [l2[ 
we experimentally confirmed that the potentials obtained 
from the simulations match the actual potentials in our 
trap with an accuracy on the 1% level. From the solver 
one can obtain normalized shape factors, Ui{x), that are 
fixed by the geometry of the trap and depend on the seg- 
ment number, i, and the distance from that segment, x. 
The total potential can then be obtained by summation 
of the voltages applied to each segment multiplied by its 
respective shape factor: 



V{i) = V 1 u l {x) + V2(t)u 2 {x) 



(1) 



ion is moved a certain distance and then moved back end- 
ing at it's initial position. V2(t) is created by a waveform 
generator and all unused segments are kept at ground. 
The duration of the transport pulse is T= 400 us, which 
is shorter than the coherence time of the qubit state 
but much longer than the inverse of the axial trap fre- 
quency, (1.35MHz) -1 . The latter condition means that 
the ion may adiabatically follow the changes to the volt- 
ages, such that the position of the ion can be obtained 
from the minimum of the total potential at any given 
time as a function of the applied segment voltages, i.e. 
2aon = Xi on (Vi, V^(i)). The total phase shift acquired by 
the ion is given by the time integral of the frequency shift 
caused by the B-ficld gradient during the transport. The 
frequency shift is proportional to the B-field, which in 
turn is proportional to the spatial position in a first or- 
der approximation. We can thus express the phase shift 
in terms of the ion position and the magnetic field gra- 
dient: 



h Sx 



[Xiony^) ^ion (0)]dt, (2) 



Here, V\ is the fixed potential on the main trapping seg- 
ment and V2(t) is a time dependent sine-shaped voltage 
that is applied to a neighboring segment, such that the 



where the Lande factor of the Si/ 2 state is given by gj = 
2. The integral in this equation is thus determined from 
the displacement of the ion, which is given by evaluating 
the minimum of the applied potential. It is then clear 
that the magnetic field gradient, SB/dx, can be extracted 
by measuring the acquired phase shift for a particular 
known displacement-time integral. 



III. RESULTS 
A. Measurement of magnetic field gradients 

The result of our experiment to measure the magnetic 
field gradient over a maximum distance of 140 urn is plot- 
ted in Fig.[2]for different currents on the Hclmholtz coils. 
Three distinct cases were measured: i) a large gradient 
corresponding to 3 A in coil 1 and 1 A in coil 2, ii) equal 
currents, 1.9 A, in each coil and iii) with currents chosen 
to minimize the gradient, yielding 1.660 A in coil 1 and 
2.122 A in coil 2. The particular values in the optimized 
last case were chosen by extrapolation from the first two 
cases, with the desire to reach a zero gradient. In all 
three cases, the offset magnetic field is close to the same 
value, i.e. around 0.7 mT, but the gradient varies over 
more than two orders of magnitude. Furthermore, the 
standard deviation of the phase measurement in Fig. [2] 
is about 2°, and from the near-perfect linear fit, we can 
conclude that the linear approximation that was made in 
Eq. ([2]) was well justified. In the figure, each phase shift 
is evaluated from a cosine fit to the fringe pattern data, 
resulting from a spin echo sequence where the phase of 
the final 7r/2-pulse is varied between and 2n in 100 
steps, as shown in the inset. The phase is then plotted 
versus the displacement-time integral (see Eq. ©), i.e. 
the representation of how long the ion stays in a magnetic 



3 



400 



350 



^300 

03 


CD 

^250 

CD 
TJ, 

^ 200 



150 
CO 
CO 
£Z 

D- 100 



50 



50, 



' VB = 3.58 • lO" 10 T/im 



which is independent of contrast and offset: 




100 200 300 

Phase (deg) 

VB = 6.74 • 10- 10 T/fim (Equal currents) 



VB = -7.9 • 10~ 12 T/fim (Minimized gradient) 



5 10 15 20 25 30 

Displacement time integral (fim ■ ms) 

FIG. 2. (color online) Recorded phase shift as a function of 
transported distance and time (see text), for three different 
coil currents. The slope gives the magnetic field gradient, 
which is listed next to each respective curve. Note how the 
optimized coil current setting has a gradient that is two or- 
ders of magnitude lower than when the coils have the same 
current. The inset shows the evaluation of a spin echo fringe 
for one particular gradient and transport setting compared to 
a non-shifted one. Each point in the inset consists of 200 ion 
interrogations. 



field that is different from its starting location, times the 
transported distance. 



B. High precision detection of small gradients 

The main application for this type of magnetic field 
inhomogeneity measurement is the creation of a low de- 
coherence environment for quantum computing, and we 
now focus on measurements of the low gradient case with 
high precision. We also aim at using as short as possible 
total measurement time, such that the method can be 
efficiently used as a calibration routine. Shortening the 
time of the acquisition can be achieved by reducing the 
number of points that are used in order to determine the 
acquired phase shift. In principle, for small phase shifts, 
it is enough to measure the spin up probability at a con- 
cluding angle of 9 = it/ 2 in the cosine curve, where the 
slope is the steepest, in order to obtain an estimation for 
the phase shift. In reality however, various mechanisms 
can cause offsets and loss of contrast of the fringe, which 
leads to the 7r/2 point getting shifted for other reasons 
than acquired phase. To accommodate for this, we ad- 
ditionally measure the baseline and the contrast of the 
fringe pattern by also probing the points 9 = and 9 = it 
respectively We can now form a normalized signal, S 1 , 



S{<j>) 



p(tt/2,0) -p(n,<j)) 
p(0, 4>) - p(n, 4>) 



(3) 



With a 5 ms detection time, one measurement of the sig- 
nal S '(</>), comprised of the three measurement points at 
<j> = {0,0.5,1}tt, thus takes 15 ms. In Eq. ©, p{6,4>) 
denotes the probability of measuring the ion in spin up 
as a function of the angle of the concluding 7r/2-pulse, 
9, and the phase shift induced by the magnetic field 
gradient, <f). In general, p(9, c/>) is described by a sinu- 
soidal function with a scaling factor and an offset, i.e. 
p{9,4>) = Ai cos(# + 4>) + A 2 . Putting this expression 
into Eq. ([3]) however, it becomes clear that the parame- 
ter S is independent of both the scaling factor and offset 
and we find that the phase shift can be obtained as 



arctan(l — 25). 



(4) 



In order to investigate the precision obtained from the 
above construction, we measured a repeated acquisition 
of these three points with and without transport, for 
small changes to the current of the coils. An example 
of this is shown in Fig. [3^, where each point is an av- 
erage of all preceding points, and each error bar is the 
deviation of the last (averaged) point with respect to the 
previous points. If the measurement is quantum shot 
noise limited, this error should converge towards a cor- 
rect estimation of the phase shift with a rate proportional 
to y/N, where N is the number of repetitions of the ex- 
periment. 
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FIG. 3. (color online) Part a) shows an example measurement 
of a small magnetic field difference of 7 nT over 19 urn (blue 
curve) in comparison to a minimized gradient (red curve), 
where each of the 100 experimental points consists of 200 ion 
interrogations. In b), an Allan type standard deviation over 
a larger data set (black dots) is plotted in comparison with 
the theoretical expectations using the values of the contrast 
going from 0.12 to 0.75 as realistic values (blue solid line) and 
using a perfect contrast (dashed red line), for the standard 
deviation coming from the quantum projection noise limit 
(see text). The y-axis scale represents the smallest magnetic 
field difference we can measure on top of the background of 
B = 0.7 mT. 
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C. Allan variance and the quantum shot noise limit 

While Fig. [3£i demonstrates how the phase shift with 
and without transport directly can be used to resolve 
small differences in the magnetic field, an Allan type eval- 
uation of the standard deviation gives a better estimate 
for the true rate of convergence. This is displayed in part 
b) of Fig. [3] In order to show that the main limitation 
of the experiment comes from quantum shot noise, we 
derive an expression for the quantum projection noise, 
caused by the projection on the spin states. This binary 
choice indicates that the uncertainty for each individual 
ion interrogation is given by the binomial distribution: 



,(0) = Vp(</») • (l-p(<l>))/N, 



(5) 



where p is the probability of measuring spin up, and N 
is the total number of individual interrogations. As each 
data point consists of three actual measurements, the 
contributions to the error from all three projections must 
be taken into account. Gaussian error propagation then 
leads to an expression for the total error of the phase 
estimation: 



0> 



*=0,7T/2,7T 



\s P (4>) 



(6) 



where the partial derivatives can be obtained from Eqs. 

and The binomial standard deviations from 

Eq. J5) are used as errors together with values for p(<j>) 
corresponding to having the three points come from a 
signal with realistic contrast going from 0.12 to 0.75 
{A\ = 0.31 and A-i = 0.44), like the one seen in the inset 
to Fig.d We then obtain a limit of A</> = 1.81 rad/V^V, 
which is plotted in Fig. \3jp, converted into a correspond- 
ing relative magnetic field difference. In addition, the fig- 
ure also displays the minimal obtainable error that could 
be reached if the spin echo contrast would be perfectly 1. 
We find that the experimental data follows the quantum 
projection noise curve, which takes the limited contrast 
into account. From the above derivation it is also clear 
that the obtained accuracy is not limited by the time it 
takes to perform the acquisition, but only by the num- 
ber of interrogations, which is substantiated by the close 
proximity to the theoretical limit. This means that tech- 
nical improvements, such as more efficient camera read 
out, can be used to reduce the time it takes for a full 
acquisition. As a reference, the plot in Fig. [3^ contains 
in total 20000 ion interrogations (for each of the three 
points), and with 5 ms detection time per interrogation 
it took only in the order of minutes to acquire. After 
about 10 4 interrogations, we reach a relative magnetic 
field sensitivity of SB/B ~ 5 • 10~ 7 . At this point the 
measured Allan deviation appears to deviate from the 
theoretically expected one. This behavior is attributed 
to long term drifts of the experimental setup arising from 
e.g. temperature variations in the magnetic field coils. 
The observed contrast loss can be explained by techni- 
cal noise from the waveform generator that applies the 



transport pulse voltage, which leads to a small amount 
of motional heating. This problem can be solved by a 
better suited multi-channel voltage source that can ad- 
dress all of the trap control segments. Time dependent 
magnetic field fluctuations faster than the spin echo time 
is further reducing the contrast, and such noise may be 
reduced by a passive shielding box of /x-metal around the 
trap setup. It can also be mentioned that it is possible 
to improve the sensitivity of the method by utilizing two 
or more ions in a spin-squeezed state, which can reduce 
the noise of the phase measurement even below the shot 
noise limit (l5j . 

Our value of the relative magnetic field sensitivity is 
comparable to other methods used for measuring niag - 
netic fields, such as squids [l6| or atomic vapors [l7l[l8j. 
Though it may be noted that the other methods target 
different applications, and while those can reach higher 
absolute sensitivity, our method relies only on a single 
ion, and can thus yield a very high spatial resolution, 
in the sub- urn range. Another potential application is 
the detection of the magnetic field from a single spin. 
Although it is currently difficult to move the ion close 
enough (~ 100 nm) to the spin for adequate absolute 
sensitivity, it is interesting to note that our method has 
a similar relative sensitivity as magnetic resonance force 
microscopy where single spin resolution has been 
reached. 



IV. CONCLUSION 

In conclusion, we have demonstrated a technique to use 
transport of a single ion to detect magnetic field gradi- 
ents with high resolution, both spatially and magnitude 
wise. Using the transport of the ion gives the natural ad- 
vantage that magnetic gradients in the entire trap, even 
regions inaccessible to the necessary lasers, can be mini- 
mized, which is expected to be very useful for ion traps 
used for scalable quantum computing. For a single trans- 
ported ion, the minimum detectable phase shift suggests 
a transported coherence time would be more than ten 
ms, which is already long enough that the gradient is 
no longer the main limitation. The expected coherence 
time is even significantly longer if we consider a logical 
qubit in a decoherence free subspace @] . This entangled 
two-ion logical state is insensitive to magnetic field inho- 
mogeneities, although a difference in the magnetic field 
between the two ions will still cause dephasing. Using 
the minimized gradient demonstrated here and assuming 
a reasonable ion separation of 4 urn, we find a coherence 
time above 1 s. This time is long compared to the time 
of qubit gates as well as the expected transport dura- 
tions, however, it is not so long that one does not realize 
the necessity of performing a compensation routine, like 
the one suggested here, for minimizing the magnetic field 
inhomogencity. 
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